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O Abstract 

<p : 

We show how parton distributions unintegrated over the parton transverse momentum, 
kt, may be generated, at NLO accuracy, from the known integrated (DGLAP-evolved) par- 
ton densities determined from global data analyses. A few numerical examples are given, 
, which demonstrate that sufficient accuracy is obtained by keeping only the LO splitting 

ly-j , functions together with the NLO integrated parton densities. However, it is important to 

keep the precise kinematics of the process, by taking the scale to be the virtuality rather 
than the transverse momentum, in order to be consistent with the calculation of the NLO 
' splitting functions. 

O ' 
> 

jx} ■ 1 Introduction 

Conventionally, hard processes involving incoming protons, such as deep-inelastic lepton-proton 
scattering, are described in terms of scale-dependent parton distribution functions (PDFs), 
a(x,/x 2 ) = xg(x,[i 2 ) or xq(x,fi 2 ). These distributions correspond to the density of partons 
in the proton with longitudinal momentum fraction x, integrated over the parton transverse 
momentum up to k t = fi. They satisfy DGLAP evolution in the factorisation scale /i, and are 
determined from global analyses of deep-inelastic and related hard-scattering data. However, 
for semi-inclusive processes, parton distributions unintegrated over k t are more appropriate. 
For example, unintegrated parton distributions play an important role in the description of 
the transverse momentum dependence of different inclusive hard processes, such as inclusive 
jet production in deep-inelastic scattering (DIS) [3], electroweak boson production [2], prompt 
photon production [3], azimuthal correlations in high-p T dijet production [I], etc. Moreover, 
the exclusive cross sections for vector meson photoproduction [5J or central exclusive diffractive 




Figure 1: A schematic diagram of inclusive jet production in DIS at LO which shows 
the approximate equality between, on the left-hand side (a), the formalism based on the 
doubly-unintegrated quark distribution, f q (x, z, k 2 , fi 2 ), where the off-shell quark has virtuality 
—k 2 /(l — z), and on the right-hand side (b), the conventional QCD approach using integrated 
parton densities, a(x,/i 2 ), where the incoming partons are on-shell. 

Higgs boson production [6] are also calculated in terms of the unintegrated parton distributions. 
In fact, so-called l k t - factorisation' was originally established [7] for heavy-quark pair production, 
so that the cross section for pp — > QQX is of the form: 

/"dxi fdx 9 fdk 2 . fdk 2 , _ „ „ 

a(pp^QQX) = J-± J-l j ^ f g ( Xl ,kl) f 9 (x 2} kl t ) a(s } M 2 } kl,kl), (1) 

where the f g are the gluon densities of the incoming protons, unintegrated over k 2 t , such that 
f g (x, k 2 )(dx / x){dk 2 / k 2 ) is the number of gluons in the longitudinal and transverse momentum 
intervals from x to x + dx and from k 2 to k 2 + dk 2 , respectively, and a is the gg — > QQ 
subprocess cross section. 

In general, the unintegrated distributions, f a (x, k 2 , /i 2 ), depend on two hard scales, k t and 
//, and so the evolution is much more complicated. The additional scale /x plays a dual role. On 
the one hand it acts as the factorisation scale, while on the other hand it controls the angular 
ordering of the partons emitted in the evolution. In Refs. [TJ |H] a prescription was given 
which allows the unintegrated distributions to be determined from the well-known integrated 
distributions. The prescription was based on the fact that due to strong k t ordering, inherent in 
DGLAP evolution, the transverse momentum of the final parton is obtained, to leading-order 
(LO) accuracy, just at the final step of the evolution. Thus the /cj-dependent distribution can be 
calculated directly from the DGLAP equation keeping only the contribution which corresponds 
to a single real emission, while all the virtual contributions from a scale equal to k t up to the 
final scale \x of the hard subprocess are resummed into a Sudakov-like T-factor. The factor T 
describes the probability that during the evolution there are no parton emissions. 

The idea is that, considering inclusive jet production in DIS, for example, the LO diagram at 
0(a em ) computed using fc f -factorisation will already include, to a good approximation, the main 
effects (which are of kinematical origin) of the conventional LO QCD diagrams at 0(a em cts) 
computed using collinear factorisation. This approximate equality is shown schematically in 
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Fig. [TJ The cross section for any hard process is then determined by convoluting the unin- 
tegrated parton distributions with the off-shell subprocess cross sections where the incoming 
partons have virtuality — k 2 . To be precise, it is necessary to also take account of the frac- 
tion z of the light-cone momentum of the parent parton carried by the 'unintegrated' parton, 
that is, to use 'doubly-unintegrated' parton distributions [T], f a (x, z, k 2 , /i 2 ), where the off- 
shell parton now has virtuality — k 2 /(l — z). The doubly (or 'fully') unintegrated distributions 
preserve the exact kinematics of the partonic subprocess (see also Ref. [H]) and we speak of 
(z, /^-factorisation. Here, we wish to extend the 'last step' LO prescription for determining 
unintegrated parton distributions to next-to-leading order (NLO)Jll 

First, in Section [2] we recall the LO prescription, then we extend it to NLO in Section [3j 
We show numerical results in Section @] and conclude in Section [51 More details of the NLO 
derivation are given in an Appendix. 

2 LO prescription for unintegrated parton distributions 

It is useful to review how LO unintegrated parton distributions, f a (x,k 2 ,n 2 ), may be calcu- 
lated from the conventional (integrated) parton densities, a(x,fi 2 ), in the case of pure DGLAP 
evolution. As usual, we adopt a physical (axial) gauge, which sums over only the transverse 
gluon polarisations, so that the ladder-type diagrams dominate the evolution. Recall that the 
number of partons in the proton with longitudinal (or, to be precise, light-cone plus^l) mo- 
mentum fraction between x and x + dx and transverse momentum k t between zero and the 
factorisation scale \x is a(x, /j, 2 )(dx/x), whereas the number of partons with longitudinal mo- 
mentum fraction between x and x + dx and transverse momentum squared between k 2 and 
k 2 + dk 2 is f a (x,k 2 , /i 2 )(dx/x)(dk 2 /k 2 ). Thus the unintegrated distributions must satisfy the 
normalisation relationn 



where a(x,fi 2 ) = xq(x,fi 2 ) or xg(x,fi 2 ). We start from the LO DGLAP equations evaluated at 
a scaled k 2 : 



1 We do not consider fully unintegrated distributions at NLO. This would require the recalculation of the 
NLO DGLAP splitting kernels in fully unintegrated form, which is a necessary ingredient for a NLO parton 
shower Monte Carlo (see Ref. [10] for work in this direction) . 

2 The plus and minus light-cone components of a parton with 4-momentum fc are fc = fc° ± fc 3 . 

3 Note that the exact value of the upper limit in Eq. ([2]), and the possible non- logarithmic tail for kt > /U, 
are beyond NLO accuracy. 

4 Usually DGLAP evolution is written in terms of the virtuality fc 2 , but at LO level this is the same. The 
difference is a NLO effect. We examine the difference in detail in Section 2] 




(2) 



dxq(x,k 2 ) = a s (k 2 ) 
d log k 2 2-k 




fdzP qb {z)b(-,k 2 )-xq{x,k 2 t ) [ d( P qq (C) , (3) 
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dxg(x,k 2 ) a s {k 



d log k 2 2n 



I dzP gb (z)b(-,k 2 )-xg(x,k 2 t ) TdC ((P gg (() + n F P qg (()) 

b=q,g Jx " J ° 



(4) 

where b(x,k 2 ) = xq(x,k 2 ) or xg(x,k 2 ) and P a b(z) are the unregulated LO DGLAP splitting 
kernels. The two terms on the right-hand sides of Eqs. ([3]) and (jlj) correspond to real emission 
and virtual contributions respectively. The virtual (loop) contributions may be resummed to 
all orders by the Sudakov form factor, 



2 k 2 2n 



T q (klfi 2 ) = exp - / -i dC P M (C) , (5) 



T,(A; 2 , /i 2 ) ee exp (- [" ^ /'dC (CP 99 (C) + ^P ?S (C))) , (6) 

y Jk 2 t K t Z7r Jo J 

which, recall, give the probability of evolving from a scale k t to a scale /i without parton 
emission. Thus, from Eqs. (J2]) (HJ) , the unintegrated distributions have the form: 



f a (x,k 2 ^ 2 ) = T a (kl^).^l Yl J dzP ab (z)b(^k 2 ). (7) 

b=q,g x 

This equation can be rewritten as 



d 

f a (x,k 2 ,fx 2 ) = — — n [a(x,k 2 )T a (k 2 ,fx 2 )] , (8) 
a log kf 

since the derivative dT a / d\ogk 2 cancels the last terms in Eqs. ([3]) and (jlj). 

This definition is meaningful for k t > fi , where /j,q ~ 1 GeV is the minimum scale for which 
DGLAP evolution of the conventional parton distributions, a(x,/j, 2 ), is valid. Integrating over 
transverse momentum up to the factorisation scale we find that 



/.(*,*?, A* 2 ) = [a(xX)T a (kl^)]l= = l 

Jul K t 



- =1*0 

'Mo * 

= a(x,fi 2 ) -a(x,fil)T a (nl,fi 2 ). (9) 



Thus, the normalisation condition Eq. ([2]) will be exactly satisfied if we assume^ 

1 



1 - ' -2 ,.2\ 



2 fa{x } k t , [1 



a(x,/4)T a (nl,ii 2 ), (10) 



so that the density of partons in the proton is constant for kt < /xo at fixed x and \i. 

So far, we have ignored the singular behaviour of the unregularised splitting kernels, P qq {z) 
and Pgg(z), at z = 1, corresponding to soft gluon emission. These soft singularities cancel 



5 A more complicated extrapolation of the contribution from k t < ^jlq, which ensures continuity of /„ at 
kt = Mo 7 w as given in Ref. [2]. 
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between the real and virtual parts of the original DGLAP equations, Eqs. ()3]) and fl3J). So 
it was enough to replace an explicit cutoff by the so-called '+' prescription. Now, when real 
emission (which changes the k t of the parton) and the virtual loop contribution (which does not 
change the kinematics of the process) are treated separately, we must introduce the infrared 
cutoff explicitly The singularities indicate a physical effect that we have not yet accounted for. 
It is the angular ordering caused by colour coherence [UJ, which implies an infrared cutoff on 
the splitting fraction z for those splitting kernels where a real gluon is emitted in the s-channel. 
Indeed, the polar angle 9 ordering 



• • • < Oi-i < Q% < Oi+i < . . . 



(11) 



implies the ordering of the ratios & = p^/pf, that is, of the rapidities rji = — (1/2) log £j = 
— logtan(#j/2). The light-cone components of the massless gluon-z momentum (pi = — ki) 
satisfy p~ = Pitlvt ■> w hile the ratio of the momenta in the proton direction is given by the 
ratio of the momentum fractions carried by the t-channel gluons (see Fig. 2 of Ref. [T]), Zi = 
kt/kt-i = x i/ x i-i- II we denote pi = Pi,t/(1 — z%) then we obtain 



&/e. 



i-l 



{Pi/ Zi-iPi-i) 2 > 1 that is Zi_i]9i_i < pi. 



(12) 



In the last step the angle is limited by the value of factorisation scale fi. If we choose // = Q 
for DIS in the Breit frame, then the last inequality reads [TJ 



(13) 



which leads to z < fi/(/JL+ k t ) and provides the inequality Qi < M . 

Thus the factorisation scale, /x, is entirely determined from the kinematics of the subprocess 
at the 'top' of the evolution ladder [TJ. So we define the infrared cutoff to be 



then the precise expressions for the unintegrated quark and gluon distributions are 



(14) 



f q (x,kt,/j, 2 ) — T q {k 2 ,jj? 



a s (k 



2vr 



dz 



X 



X 



X 



P qq (z) -q - k 2 t 9(1 - A - z) + P qg (z) -g -, k 



x 



and 



2tt 



dz 



(15) 

P 9q (z) % (f , k\) + P gg (z) X -g (J kt) 9(1 - A - 

(16) 

By unitarity the same form of the cutoff, A(k 4 ), must be chosen in the virtual part. Thus 
we insert 9(1 — A — Q into the Sudakov factor for those splitting functions where a gluon is 



ft i, 2 2\ rp /<,2 2\ a s(kt) 

Jg{X, K t , \l ) — lg{K t , fi 
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emitted in the s-channel and 0(£ — A) where a gluon is emitted in the t-channeljj Then 

T q (kl^) = exp (- ^1 jfV^COeCl - A - C )j . (17) 

Recall that the exponent of the gluon Sudakov factor was already simplified by exploiting the 
symmetry P qg (l — C) = Pqg((), so that the gluon Sudakov factor is 

T g {kl /i 2 ) = exp {- jT ^ ^1 j\( [ C P OT (C)e(l - A - C)6(C - A) + n F P qg {Q ] J , 

(18) 

where np is the active number of quark-antiquark flavours into which the gluon may split. 

Note that in the expression for the unintegrated distribution, Eq. (jSJ), the derivative of a 
with respect to log k\ gives the right-hand side of the usual DGLAP equation, which contains 
both the real emission (which we are looking for) and the virtual loop contribution. However 
the virtual contribution is cancelled by the derivative of T a with respect to log k\ . 

In the next section we introduce an analogous prescription for unintegrated parton distri- 
butions which may be justified at NLO accuracy. 



3 Unintegrated parton distributions at NLO 

First, we discuss the structure of the NLO contribution. The original definition was that the 
amplitude was decomposed into the perturbative sum 



oo oo 



A = ^^A n , fc 4KlogQ 2 )", (19) 

k n=l 

where k — 0, 1, 2, . . . denote the LO, NLO, NNLO, . . . contributions. However, nowadays the 
analogous decomposition is used for the anomalous dimensions, 7^, which describe the evolution 
in terms of the OPE operators, e.g. Mellin moments, M^, 

M N (Q) = M N (Q )(Q 2 /Qly», (20) 

rather than for the amplitude. That is, Mjv is written in the form 



M N (Q) = M N (Q ) exp ( log(Q 2 /Q 2 ) £ c N k a k s ) = M N (Q ) £ ( log(Q 2 /Q 2 ) £ 1c 



n 

..V . . I. 



k a s 



(21) 

Thus the Feynman diagrams which describe DGLAP evolution have the structure shown in 



6 The lower cutoff is beyond the LO accuracy and was simply introduced to make the formulation more 
symmetric. It is not included in the NLO prescription of Section [3] 
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• p(0) 
ki+l - ' ~ - . 

<;><O p(1) 

ki " - — - 



Figure 2: The ladder diagram describing DGLAP evolution at NLO. is the appropriate 
NLO coefficient function, and P^ and P^ are the appropriate LO and NLO splitting functions. 



Fig. |2j We have strong ordering of the transverse momenta, k^ 3> fcj-i,t, between the splitting 
functions, but no ordering in the transverse momenta of the t-channel partons which enter the 
NLO (NNLO, . . .) splitting functions. 

Here, we discuss NLO. If we denote the transverse momentum in the NLO loop as k' it , then 
we have the ordering 

ht < K,t ~ h+i,t < h+2,t- (22) 

The 'parton' of the /c t -factorisation approach at NLO is the t-channel quark or gluon placed 
between the coefficient and splitting functions. That is, the parton labelled as k n in Fig. |2j 
There should be strong k t ordering between the loops which correspond to the NLO coefficient 
function and to the NLO splitting function. Otherwise the uppermost two loops should be 
assigned to the NNLO coefficient function. 

We seek the prescription which gives the unintegrated k n>t distribution of the parton to 
NLO accuracy. To do this we have to extend the formalism of the last section in two steps. 
First, we have to replace the LO splitting functions P a b{z) by the corresponding 'LO+NLO' 
splitting functions [121 13] an d the LO 'global' parton densities a by the NLO parton densities. 
Secondly, we must take care of the precise value of the scale at which the parton densities a are 
measured, and of the limits in the integrations over the terms which include the LO splitting 
functions. In the DGLAP evolution equation the current scale k 2 is not exactly equal to the 
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parton transverse momentum k t , instead we have^ 

k 2 

k 2 = (23) 

At LO level, or at very small z, this scale difference is negligible. However, to reach NLO 
accuracy we have to account for this effect, at least in the LO part of the splitting functions. 
Thus, recalling Eq. ((7j), at NLO we have 

f™°(x,k>y) = /\i*T B (jfcV) y: ^ o+i) w ^ nlo {- z ,# = S^) Q ( i -^v), 

X b=q,g ^ ' 

(24) 

where P(° +1 ) = P<® + (a s /2n)P ( - 1 \ and 

= ^°(«) - ©(* " (1 " A)) F« <U Pa6 (s), (25) 



with strength of the z — > 1 singularity in given by [13] 

= F« = -C F (T R N F f+C A (%-%)), (26) 

= 2C A , F$ = -2C A (T R N F f + C A (f-%)). (27) 

Here % — 0, 1 denote the LO and NLO contributions, respectively, and p qq (z) = (1 + -2 2 )/(l — 2) 
and Pgg(z) — z/(l — z) + (1 — z)/z + z(l — z). More details are given in the Appendix. 

The last term in Eq. (1231) accounts for the coherence and eliminates the 1/(1 — z) singularity 
in the splitting function caused by the emission of one soft gluon, which violates angular ordering 
Eqs. (I13II14I) . Recall that at NLO, the only strong singularity in P a b comes from the vertex and 
self-energy corrections to single gluon emission, see Table 1 of Ref . [12] . The contribution from 
the emission of two 'real' gluons has (after the usual subtraction of that generated by two LO 
splitting functions) no 1/(1 — z) singularity, but only the 'soft' integrable term proportional to 
log(l — z). Note that, in contrast to Eq. (j7|), at NLO T a and a$ occur inside the z integration 
as k 2 = k 2 /(l-z). 

The final point is to consider the Sudakov factors^] T a which resum the virtual DGLAP 
contributions during the evolution from k 2 to /1 2 . At the NLO we get 

T q (k 2 ,^ 2 ) =exp Cft° +1) (C)+Pr i} (C)]) , (28) 



Eq. (|23|) assumes that the mass of the partonic system produced by the splitting is zero. This is true at 
LO, but not at NLO, where a pair of massless partons may be created. However, to account for non-zero mass 
of the pair, in the NLO splitting function, is a NNLO correction, which is beyond our NLO accuracy. 

8 Note that, at NLO, the splitting function P qq contains terms which change the flavour of the quark, P qiqj ■ 
So now T q of Eqs. (|17I28I) must include a sum over the flavour of the new quark i, and similarly for the sum 
over the initial quark b — qj in Eq. (l24)l . 



8 



T g {k 2 ^ 2 ) = exp (-[^^ J\ m^iO +2^ g 0+1) (0]) , (29) 
where, at NLO, we must take p(° +1 ) = pM> + (a s /2ir)P^ and 

A = — — with Kt = Jk 2 {1-0- (30) 

fJL+Kt 

The theta functions in Eq. (|25|) again provide the correct angular ordering for soft gluon emis- 
sion. If n t <C /i, the effect of A in the non-singular part of P is negligible. The region of K t close 
to fj, should be specified by the convention of how to separate the NLO coefficient function (that 
is, the hard matrix element) and the unintegrated parton distributions of the fc t -factorisation 
approach. In particular, the diagram corresponding to the parton self-energy, or to large an- 
gle gluon emission, could be assigned either way; that is, to the parton distribution or to the 
coefficient function. Our O functions correspond to strong angular ordering. Our prescription 
means that the contribution coming from angles 6 > 9^ should be included in the coefficient 
function, while the lower 6 domain should be assigned to the unintegrated parton distribution. 
The angle 6^ is defined below Eq. ( TT3l) . 

4 Numerical results at NLO 

For simplicity, we concentrate on the singlet evolution with np = 3 and we use the corre- 
sponding MSTW 2008 PDFs [14] and as evolved with three fixed flavours. First, in Fig. [3] 
we show the DGLAP splitting functions at LO and NLO given by Eq. fl25|) . It is seen that 
the NLO corrections to the splitting functions are relatively small in comparison with the LO 
contributions. The largest NLO correction is perhaps in P qg at small z, where at LO there is 
no 1/z term as z — > 0. In Fig. [4] we present the unintegrated PDFs as a function of k 2 at 
/i 2 = 10 4 GeV 2 for x = 0.1, 0.01, 0.001 and 0.0001. The corresponding plot at fx 2 = 100 GeV 2 
is shown in Fig. [51 Here, q is the quark singlet distribution, q = u + d + s + u + d + s. Since 
the LO and NLO splitting functions are similar, the unintegrated PDFs calculated using the 
full NLO framework are close to the results keeping only the LO part of the splitting functions: 
compare the dotted and dashed curves in Figs. 0] and [51 If we were to use the LO (rather than 
NLO) integrated PDFs, then we obtain the unintegrated distributions shown by the continuous 
curves. We see that there is a sizeable enhancement of the unintegrated gluon at very small 
x, that is x < 10~ 3 . This simply reflects the well-known difference between the integrated LO 
and NLO gluon distributions at small x0 The kink in the k t distributions at relatively large 
k t is due to the presence of two cutoffs in the expressions for the unintegrated distributions. 
One cutoff, 0(1 — A — z), accounts for the coherence effect, which leads to angular ordering. 
The other cutoff is due to the bound on the virtuality, k 2 < fi 2 , in the DGLAP evolution. For 
k t <C /i, the angular ordering of the soft gluon emissions is the stronger bound, while at large k t 
the bound on the virtuality takes over, with the transition point at kt = (1/2) (v^— 1) ^ — 0.6/x. 

9 The LO integrated gluon distributions obtained in the global analyses are larger in order to compensate for 
the absence of the l/z pole in P qg at LO and the absence of a photon-gluon coefficient function at LO. 
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Figure 3: DGLAP splitting functions, zP a b(z), given by Eq. ( 1251) . at LO and NLO, after the 
subtraction for z > ^/(fi + k t ) due to angular ordering, where we take fi 2 = 10 4 GeV 2 and 
k 2 t = 100 GeV 2 . 
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Unintegrated PDFs at \i = 10 4 GeV 2 
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Figure 4: Unintegrated parton distributions, f a (x, k 2 ,/i 2 ), given by Eq. (124|) . at ji 2 = 10 4 GeV 2 , 
for different orders of integrated PDF and 'last-step' splitting function. We use MSTW 2008 [14) 
integrated PDFs to generate the numerical predictions throughout this paper, except for Fig. [7] 
which shows the relative insensitivity to the choice of input PDF set. 
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Unintegrated PDFs at \i = 10 2 GeV 2 
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Figure 5: Unintegrated parton distributions, f a (x, fc 2 ,/z 2 ), given by Eq. (GIJ, at /i 2 = 100 GeV 2 , 
for different orders of integrated PDF and 'last-step' splitting function. 
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Unintegrated PDFs at \i = 10 4 GeV 



NLO PDF, LO P , 

ab 

— Same with scale k. 
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Figure 6: Unintegrated parton distributions, f a (x, k 2 ,/i 2 ), given by Eq. (124|) . at ji 2 = 10 4 GeV 2 , 
compared to a prescription where the scale is taken as parton transverse momentum k 2 rather 
than virtuality k 2 = k 2 /(l — z). We use NLO PDFs with LO P a b in the last step in both cases. 
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As noted above, the results obtained using the full NLO prescription (the dashed curves in 
Figs. |H and ED are very close to that obtained simply using the LO splitting functions (dotted 
curves). However, they differ from the results obtained from the old LO prescription [TJ [8], 
described in Section 2, even when the NLO integrated PDFs are used in both cases. The 
comparison is shown in Fig. [6j The difference arises because in the old LO prescription we 
did not care about the precise scale. At LO accuracy, scales k 2 and k 2 are both acceptable. 
However, at NLO accuracy, we must impose the correct scale, namely k 2 = k 2 /(l — z), as in 
Eq. ( )23|) . The main effect is the larger values of the Sudakov factors, T, due to larger lower 
limits in the integrals in Eqs. ( )28|29|) leading to a smaller power in the exponent. This enlarges 
the unintegrated PDFs for k t <C /i. For large k t , however, the role of the T factor is less 
important: T is already close to 1. On the other hand, a larger scale in the integrated PDFs 
and in as leads to somewhat smaller unintegrated densities, especially at relatively large x. 
This is seen in Fig. |6l where the dashed curve corresponds to the results obtained using the 
scale k\ , instead of k 2 . 

We compare results obtained using MSTW 2008 [H] and CTEQ6.6 [15] PDFs in Fig. 
which shows that the predictions for NLO unintegrated PDFs are quite insensitive to different 
choices of the input integrated PDF set. 

Usually the unintegrated distributions are defined as the density of partons at fixed k t , 
and fc r factorisation is used to calculate the observables. In some specific cases, it may be 
advantageous to fix the virtuality k 2 rather than k 2 t of the parton0 For completeness, we show 
an example of such unintegrated distributions in Fig. [SJ Here we have used the LO splitting 
functions in the 'last step' and NLO integrated PDFs, so the curves should be compared with 
the dotted lines in Fig. |H Note that if we were to work in terms of k 2 then we would not 
have the last ©-function, 0(1 — z — k 2 /fi 2 ), which provides the inequality k 2 < fi 2 in Eq. (|24l) . 
Instead we have to introduce an infrared cutoff k t > k = 1 GeV, that is 0(A; 2 (1 — z) — k%). 
(Physically, at large distances the system becomes colourless and so the kt is limited by the 
uncertainty principle.) This infrared cutoff decreases the parton densities with fixed virtuality 
at low k 2 , whereas the usual angular-ordering ©-function takes over for large k 2 , again leading 
to a kink at the transition point where k 2 = k\ + k \i. 



5 Conclusions 

We have presented a prescription to generate unintegrated PDFs with NLO accuracy based 
on NLO integrated PDFs known from global parton analyses. The prescription comes from a 
detailed consideration of the last step of the DGLAP evolution and accounts for the precise 

10 That is, in Eqs. (jTJ) and ([2]), fc 2 is replaced by fc 2 , and we may speak of 'fc-factorisation' rather than 
fct-factorisation. Of course, in this case new coefficient functions would have to be defined and calculated 
corresponding to the unintegrated parton distributions at fixed fc 2 . Note that 'factorisation' is just a convenient 
way to write some amplitude or cross section in terms of the 'Feynman diagram blocks' (in the physical planar 
gauge); see, for example, Fig. [5] 
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Unintegrated PDFs at \x 2 = 10 4 GeV 2 
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Figure 7: Unintegrated parton distributions, f a (x, fc 2 ,/i 2 ), given by Eq. (1241) . at /i 2 = 10 4 GeV 2 
for different input NLO PDFs: MSTW [H] and CTEQ [15] (both using LO P ab in the last step). 
For comparison to CTEQ6.6, here we use the standard MSTW 2008 PDFs and a$ evolved with 
rip = 5, leading to a slight difference compared to the other plots in this paper where the PDFs 
and as are instead evolved with rip = 3. 
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Figure 8: Unintegrated parton distributions, f a (x,k 2 , /i 2 ), given by Eq. ( lA.lljl . at fi 2 = 10 4 
GeV 2 , defined as a function of the virtuality k 2 = k\j (1 — z). We use NLO PDFs with LO P a b 
in the last step. 



16 



kinematics and the coherence effect (angular ordering) of the soft gluon emissions. We compared 
the NLO and LO prescriptions and demonstrated that, to good accuracy, the NLO prescription 
can be simplified to keep only the LO splitting functions. Nevertheless, the difference between 
the results obtained using the new NLO and the old LO prescriptions is not negligible (see 
Fig. [6]). It turns out that it is important to account for the more precise kinematics in the 
LO splitting functions by taking the scale to be the virtuality k 2 rather than the transverse 
momentum k 2 . Since the virtuality k 2 > k 2 , the Sudakov suppression, which is driven at NLO 
by the value of k 2 , is weaker now. This leads to larger unintegrated PDFs in the region k 2 <C fi 2 . 

Finally, a comment on the implications of these new predictions. As a topical example 
we consider the use of the new unintegrated gluons in the prediction of the cross section of 
exclusive Higgs production, pp — » p + H + p, at the LHC. The conclusion at the end of the last 
paragraph indicates that the exclusive Higgs cross section calculated to NLO accuracy should 
be larger than that calculated using the old LO prescription. Note that the main contribution 
to the cross section [5] comes from the region k 2 ~ 4 GeV 2 < M 2 H and x ~ 0.01, where the 
NLO prescription enhances the gluon density. Since the predicted cross section is proportional 
to fg, this enhancement has important implications for the viability of the exclusive Higgs 
signal at the LHC. This is not the final effect, since diffractive Higgs production is driven 
by generalised skewed (not diagonal) unintegrated gluons. The NLO skewed unintegrated 
PDFs can, in principle, be obtained by applying an analogous procedure to the NLO evolution 
of generalised PDFs. At NLO we anticipate that the generalised PDFs will have a similar 
kinematical enhancement, but not so strong for the exclusive Higgs cross section, since the 
skewed distribution, f g {x,x' , k 2 , /i 2 ), in the relevant region, x' <C x, is proportional to y/T and 
not to the full Sudakov factor T. 



Appendix 

Here we describe in more detail how the NLO prescription Eq. fl24|) follows from NLO DGLAP 
evolution. The DGLAP splitting functions for singlet evolution up to NLO can be written as 

pr , w=^ ) w+^^ , w- (a.i) 

271 

The splitting function at each order (i = 0, 1) can be written as an unregularised part and a 
part proportional to 5(1 — z), i.e. 

r$(z) = P!$(z)-5{l-z)K®6 ab , (A.2) 

where the soft singularities as z — y 1 cancel between the two terms after convoluting with 
the parton distributions. The coefficients of 5(1 — z) can be obtained from conservation of 
momentum fraction, i.e. 

E Ac cn ( : } (o = o =► k® = e Ac ftoCO- ( A - 3 ) 
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We define the Sudakov form factor as 



T a (k ,y. ) = exp - / — — — ^ / dC C^a (0 



(A.4) 



where 



^"w^w (a.5) 

The DGLAP equation for the integrated PDFs, a(x, k 2 ), can be written as 

[z)b(-,k 2 ) (A.6) 



da{x,k 2 ) _ a s {k 2 ) ^ f 1 , .(o+i) f l , p 
aiogfc 2 ~ 2vr 2^ ./, d2/ab [Z) b \z 



b=q,g 

a s (k 2 ^ ' < ] 



2tt 

b=q,g 



Z([ 1 dzP^ 1 Hz)b(*h>)-a(x,k*) fdCC^riC)) (A.7) 
_asijf)_ f 1 (o+i), s ,fx 2 \ a(x,k 2 ) dT a (k 2 ,fi 2 ) 

~ £j x dzP °* [z) Hi-* ; ~ T a (k 2 ,n 2 ) dhgk 2 ■ (A - 8) 

The unintegrated parton distribution, as a function of the light-cone momentum fraction x, 
the virtuality k 2 and the factorisation scale /i 2 , is defined for k 2 < fi 2 as 

Mx,k\p. 2 ) = [a(x,k 2 ) T„(k\p. 2 )] (A.9) 

To get the unintegrated distribution as a function of transverse momentum we need to move 
the Sudakov factor and as inside the integral and account for the scale dependence Eq. ( !23l) . 
k 2 = k 2 /(l — z) using LO kinematics. Then we obtain Eq. f[2~4"|) . i.e. 



fa{ x j k 2 , /i 2 ^ 



[\zT a (k 2 ^ 2 )^J2 P ^ 1} ^ b & k2 ) ®(» 2 -k 2 ). (A.12) 

Jx ^ b=q,g 



Accounting for angular ordering, which additionally constrains the kinematics, in Eqs. ( 1A.4I) . 
f lA~TT|) and (1AT2]) we replace P by P, as given in Eq. 
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